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Abstract 

We consider dynamical systems that are described by fractional power of coordinates and mo- 
menta. The fractional powers can be considered as a convenient way to describe systems in the 
fractional dimension space. For the usual space the fractional systems are non-Hamiltonian. Gener- 
alized transport equation is derived from Liouville and Bogoliubov equations for fractional systems. 
Fractional generalization of average values and reduced distribution functions are defined. Hydro- 
dynamic equations for fractional systems are derived from the generalized transport equation. 
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1 Introduction 

Integrals and derivatives of fractional order [1, 2] have found many applications in recent studies in 
physics [3, 4, 5]. The interest in fractional equations has been growing continually during the last few 
years because of numerous applications: kinetic theories of systems with chaotic and pseudochaotic 
dynamics [6, 7, 8, 9]; dynamics in a complex or porous media [10, 11, 12, 13, 14, 15]; electrodynamics 
[19, 20], dynamical systems [21, 22, 23, 24, 25, 26] and many others. 

The new type of problems has rapidly increased interest in areas in which the fractal features 
of a process or the medium impose a necessity of applying the tools that are non-traditional in 
"regular" smooth physical equations. In many problems the real fractal structure of matter can be 
disregarded and the medium can be replaced by some smoothed continuous mathematical model where 
fractional integrals appear [11, 12, 14, 20]. The order of fractional integral is equal to the fractal mass 
dimension of medium and in this way one can take into account the fractality of media. It was proved 
that integrals on net of fractals can be approximated by fractional integrals [16] and that fractional 
integrals can be considered as integrals over the space with fractional dimension up to a numerical 
factor [17, 18]. 

It is known that Bogoliubov equations can be derived from the Liouville equation and the definition 
of average value Refs. [29, 30, 31, 32]. In Refs. [17], the Liouville equation for fractional systems 
is derived from the fractional normalization condition. In [18], the first Bogoliubov equation for 
fractional systems is derived from the Liouville equation. It is known that hydrodynamic equations 
can be derived from generalized transport equation [28] which is directly derived from first Bogoliubov 
equation [29, 31, 33]. In this paper, the hydrodynamic equations for fractional systems are derived 
from the generalized Enskog transport equations. 

In Sec. 2, the fractional average values and some notations are considered. In Sec. 3, we define the 
reduced one-particle and two-particle distribution functions. In Sec. 4, the Liouville and first fractional 



Bogoliubov equations for fractional systems are considered. In Sec. 5, we derive the fractional analog 
of the Enskog transport equation. In Sec. 6, the hydrodynamic equations for fractional systems are 
derived from the generalized transport equations. Finally, a short conclusion is given in Sec. 7. 

2 Fractional Systems 

Let us consider a dynamical system that is described by the coordinates qk and momenta pk that 
satisfy the equations of motion: 

^ = S' ^ = ^^(^"'^"'*"^ k = l,...,n. (1) 

Let us introduce the dimensionless variables 

Qk Pk _, t fk 

Qk = — , Pk = — , * = IT' ^k = -^, 
go Po to Fq 

where qq is a characteristic scale in the configuration space; pq is a characteristic momentum, Fq is a 
characteristic value of force, and to is a typical time. Then Eq. (1) has the form 

dqk Pk dqk , . , . 

-Z- = — > -^ = ZFk{q,p,t), (2) 

at m at 

where 

m = Mqo/topo, Z = toFo/pQ. (3) 

are dimensionless parameters. Using the dimensionless variables {q,p,t), we can generalize Eq. (2) 
for arbitrary powers of qk and pk'- 

^ = ^, ^ = ZFk{q-,p-,t), (4) 

at m at 

where 

qk=P{Q){Qkr = sgn{qk)\qkr, (5) 

pt = miPkr = sgn{pk)\pkr. (6) 

Here k = 1, ...,n, and (3{x) = [sgn{x)]°^~^ . The function sgn{x) is equal to +1 for x > 0, and —1 for 
X < 0. 

A system is called a fractional system if the phase space is described by the fractional powers of 
coordinates (5) and momenta (6). We can describe the fractional systems in the usual phase space 
{q.,p) and in the fractional phase space ((?",p"). In the second case, the equations of motion for the 
fractional systems are more simple. Therefore we use the fractional phase space. The fractional space 
is considered as a space with the fractional measure that is used in the fractional integrals. 

The generalization of Hamiltonian system is described by 



dql_ dH_ dp% _ dH 
~df~dp^' 'dr~~dq^' 

where H is an analog of the Hamiltonian. Using the fractional power Poisson brackets 



(7) 



_^ dA_dB__dA_dB\ 



we rewrite Eq. (7) in the from 



'^ = {ql,H}^, ^ = {p-,HU. (9) 



These equations describe the system in the fractional phase space {q°',p°')- For the usual phase space 
{q,p), Eq. (7) has the form 



dqk (qkPk) " dH dpk {qkPk) " dH 



(10) 



dt o? dpk ' dt a^ dqk ' 

and these systems are non-Hamiltonian systems. A classical system is called Hamiltonian if the right- 
hand sides of the equations 

-^ = 9k{q,p), -^ = fk{q,p) (11) 

satisfy the Helmholtz conditions [27]: 

dgk _ dgi _ dgk _ dfi _ dfk _ dfi _ , . 

dpi dpk'^' dqi dpk'^' dqi dqk'^- ^^'' 

It is easy to prove these conditions are not satisfied for Eq. (10). Therefore the dymanical system 
(10) is a non-Hamiltonian system. The fractional phase space allows us to write Eq. (9) in the simple 
form (7) and describe some non-Hamiltonian systems as a Hamiltonian systems in generalized space. 
If dq'^/dt = p^,/m, then an analog of Hamiltonian 

The omega function for system (11) is defined by 

^-^^oqk dpk^ 

and describes the velocity of phase volume change. If fi < 0, then the system is called a dissipative 
system. If Q 7^ 0, then the system is a generalized dissipative system. For system (10), the omega 
function (14) is not equal to zero, and the systems is the general dissipative system. 

It is not hard to prove that Hamiltonian (13) is connected with the non-Gaussian statistics. Dis- 
sipative and non-Hamiltonian systems can have the canonical Gibbs distribution as a solution of the 
stationary Liouville equations [40] . Using [40] , it is easy to prove that some of fractional systems can 
have fractional Gibbs distribution 

p{q,p) = exp [T - Ha,/3{q,p)]/kT, (15) 

as a solution of Liouville equation for fractional systems [17]. 

3 Fractional Average Values and Reduced Distributions 

3.1 Fractional average values for configuration space 

Let us derive the fractional generalization of average value of classical observable A(q,p). For config- 
uration space, the usual average value is 

/+00 
A{x)p{x)dx, (16) 

-00 



and can be written 

A{x)p{x)dx + / A{x)p{x)dx. (17) 

-oo J y 

Using 

^-^^^^^ T{a)jy (x-y)i-"' 
we rewtite Eq. (17) in the form 

<A>^={l\Ap){y) + {llAp){y). 

The fractional generahzation of this equation is 

<A>^={IlAp){y) + {I'lAp){y). (18) 

We can rewrite Eq. (18) in the form 

1 f°^ 
<A>^=- [{Ap){y-x) + {Ap){y + x)]dpa{x), (19) 

•^ J —oo 

where 

(i^a(x) = = x" = s5n(x)|x|°. (20) 

1 (a) al (aj 

Equation (18) defines the fractional generalization of the average value for coordinate space. 

3.2 Fractional average values for phase space 

Let us introduce some notations to define the fractional average value for phase space. Tilde operators 

J{...,Xk,...) = -(/(..., Xfc-Xfc,...) +/(...,Xfe + Xfc,...)) 



J-XkJ 



allows us to rewrite 



- (^A{q' -q,p' - p, t)p{q' - q,p' - p, t) + A[q' + q,p' - p, t)p{q' + q,p' - p, t) + 

+A{q' -q,p' + p, t)p{q' - q,p' + p,t) + A{q' + q,p + p, t)p{q' + q,p' +p,tyj 

in the simple form 

TqTp{A{q,p,t)p{q,p,t)). 

For k particle with coordinates qks and momenta pks, where s = 1, ...,m, we define the operator 

For the n-particle system phase space, we use 

T[l,...,n] =r[l]. ..T[n]. 
Let us define the integral operators /"^ by 

r+co 

I^Jixk) = / f{xk)dpUxk), (21) 



then Eq. (19) has the form 

<A>^=i^T^A{x)p{x). 

For A;-particle we use the operator 

faih.] _ fa fa fa fa 

^ ' qkl Pkl" qkm Pkml 

such that 

-^"W/(qfc,Pfc) = / /(qfc,Pfc)t?^a(qfc,Pfe), (22) 

where ducti'^ki'Pk) is an elementary 2r7T.-dimensional phase volume 

(i^Q(qfc, Pfc) = (ar(a))"^™d^^i A dpli A ... A dq^n, A dpl^n- 
For the n-particle system, we use 

/°[l,...,n] =/"[l].../°[n]. 

Using the suggested notations, we can define the fractional generalization of the average value for 
n-particle by 

<A>„=/°[l,...,n]T[l,...,n]^/?„. (23) 

In the simple case (n = m = 1), we have 

OO /"CO 



<A>a= / dfi^{q,p)TgTpA{q,p)p{q,p). (24) 

The fractional generalization of normalization condition [17] can be written by 

< 1 >a= 1. 

3.3 Reduced distribution functions 

Let us consider a classical system with fixed number n of identical particles. Suppose k particle 
is described by the dimensionless generalized coordinates q^s and generalized momenta pks, where 
s = l,...,m. We use the notations qfc = (qki, ■■■,Qkm) and p^ = {pki, ■■■,Pkm)- The state of this system 
is described by dimensionless n-particle distribution function />„ in the 2mra-dimensional phase space 

/>„(q,p,t) = p(qi,pi,...,q„,p„,t). (25) 

We assume that function (25) is invariant under the permutations of identical particles [34]: 

p(-",qfc,Pfc,---,q/,p/, ••-,*) = p{--,(ii,pi,-,<ik,Pk,-,t). 

Then the average values can be simplified [34]. Using the tilde distribution functions 

Pn{<i, P, t) = r[l, ..., n]/)„(q, p, t), (26) 

we define 

Pi(q,P,*) =p(qi,Pi,i) = I"['^,-,n]pn{<i,P,t), (27) 

which is one-particle reduced distribution function. Obviously, that pi satisfies the normalization 
condition [17]: 

/"[l]/5i(q,p,t) = l. (28) 

Two-particle reduced distribution function p2 is defined by the fractional integration of /5„ over all q^ 
and pfc, except k = 1,2: 

^2(q,P,i) =p(qi,Pi,q2,P2,t) = /"[3, ...,n]p„(q,p,t). (29) 



4 Liouville and Bogoliubov Equations for Fractional Systems 

Let us consider the Hamilton's equations for n-particle system in the form 

^ = G',iq^,pn, ^=ZF^iq'',p",t), (30) 

where Z is defined in Eq. (3). The evolution of /?„ is described by the Liouville equation [17] for 
fractional system 

^ + n^pn = 0. (31) 

This equation can be derived [17] from the fractional normalization condition 

r[l,...,n]p„(q,p,t) = l. (32) 

In Eq. (31) the derivative d/dt is a total time derivative 

, o n,m , /-, n,m , ^ 

d c* , Y^ ggfcs o Y^ dpks o 



dt dt ^^ dt dqks ^^^ dt dpks 
that can be written for the fractional powers 



The a-omega function is 



where 



, rj n,m o n,inn n 



""= i;({Gj,pJ.}« + Z{&,F'}<.). (3-1) 

k,s=l 

_ V^ ( dA^dB_ _ dA_dB\ 



Using Eqs. (34) and (33), we get Eq. (31) in the form 

dpr. 



dt 
where A„ is Liouville operator: 



AnPn, (36) 



The Bogoliubov equations [29, 30, 31, 32] describe the evolution of the reduced distribution func- 
tions, and can be derived from the Liouville equation. In Ref. [18], we derive the first fractional 
Bogoliubov equation from Eq. (36) : 

^ " d{Glh) ^^f^ djFj^p,) ^ 
Here, /(/O2) is a term with two-particle reduced distribution function, 

iih) = -T.^"[^\Pl"h. (39) 



Equation (38) is called a first Bogoliubov equation for fractional systems. 

The physical meaning of the term I{p2) is following: The term I(p2)dfj,aic[,p) is a velocity of 
particle number change in 4m-dimensional elementary phase volume dna{<ii, P2 5 q2) P2)- This change 
is caused by the interactions between particles. If a = 1, then we have the first Bogoliubov equation 
for non-Hamiltonian systems. 

5 Transport Equation for Fractional Systems 

Let us define the coordinate distribution (the density of number of particles) by the equation 

n(q,t)=r[p]pi(q,p,t), (40) 

where /[p] is a fractional integration over the momenta 

m 

^"[p] = n4"=4"-c- (41) 

s=l 

We can define the local mean values by 

(^)p,„ = (A)p,„(q,t) = ^_/"[p]yl(q,p)pi(q,p,t). (42) 

n(q,t) 

In the general case, 

{A)p,a /< A >„, (43) 

and (l)p,a = 1. The fractional average value < A >q, is connected with the mean value (A)p,a by 

<.4>«=/"[q]n(q,t)(^)p,a, 

where /[q] is a fractional integration over the coordinates q = (qi, ...,q„i). 
Fractional analog of mean local velocity is 

Vs{q,t) = — i-r[p]G,(q,p)pi(q,p,t), (44) 

n(q,t) 

i.e. 

Vs{q,t) = {Gs)p,a, (45) 

where Gs = (^^(q, p) is defined by (30). We can consider Gg as a fractional generalization of the 
velocity that has the form 

Gs = $7- (46) 

For the fractional generalization of kinetic energy of relative motion 

Y.^{p''s-Vs)\ (47) 



.=1 2 



we define density of this energy by 

M 



i?(q,t) = :^r[p](G-V)Vi(q,P,t), (48) 



where G = p"/m, and V = V(q, t). The local temperature r(q, t) is defined by the mean kinetic 
energy of relative motion: 

T{ci,t)= 7 ' 

3kBn[(i,t) 

To derive Enskog transport equation for fractional systems, we multiply both sides of (38) by the 
observable ^(p), and integrate with respect to momenta. The first and second terms of left hand side 
of Eq. (38) are transformed by 

i[p]A^ = ^^1p]M = |n(q,t)(A),,„(q,t), (49) 

^[P]^^^ = ^^^P^^G^P^ = ^nicMAGs),,^. (50) 

Integrating by part the third term of Eq. (38) and using the boundary condition 

lirn pi(q,p,t) =0, (51) 

we get 

\+oo . dA dA 

AFsPi) _^- mFsh^ = -n(q,t)(F, — ),,„. (52) 

Then we use the usual assumption 

/[p]A(p)/(p2) = 0. (53) 

ioi A = M, A = p'^ and A = p^". 

Finally, we obtain the Enskog transport equation for fractional systems: 

d d dA 

- (n(q,t)(A)p,„) + — (n(q,t)(AG)p,«) = n(q, t)Z(F^)p,„, (54) 

where 

6 Hydrodynamic Equation for Fractional Systems 

Let us consider the special cases of transport equation (54) for 

A = M, A=p<^ = MGs, A=^. 

If we use A = M, then (54) gives 

d d 

^/5M(q, t) + ^PM(q, t)(Gs)p,a(q, t) = 0, (55) 

where pM is mass density 

PAf(q,t)=Mn(q,t). (56) 



For A = p^ = MGs , we obtain 



d d 

-Q^PMVl{q^,t) + -g-^PM{GsGi)p^a = PM{(l,t)Z{Fi)p^a, (57) 



where we use Eqs. (46), (56), and the relation 

dp<^ 

For A = p2"/2M, we get 



M5is. 



d p^" d 1 

— PM(q,i)(^^)p,a + -Q-^pM{^,t){-GfGs)p,a = PM{^,t)Z {FsGs)p,a- (58) 

Here {FsGs)p,a is a local mean value. 

Let us define the deviation of velocity from its mean value by 

a(q, t) = Gs- Vsin, t) = Gs- {Gs)p,a. (59) 

Substituting Gs = Vs + Gg in the kinetic energy tensor {GsGi)p^a, we get 

{GsGi)p,a = VsVs - {GsGi)p,a, (60) 

where we use {Cs)p,a = 0. From Eq. (60), we have 

.p2° MV2 MC2 ^ ^ 

{^)p,a = ^- + -^- (61) 

The tensor of internal stress 

Psl = {CsCl)p^a, (62) 

can be represented as the sum 

Ps/(q,t) = 5siP{q,t) +TTsi{q,t), 7rss(q, t) = 0, 
where iTsi{(i,t) is the tensor of viscous stress. Then 

^{GfGs)p,a = ^{{Vi + GifiVs + Gs))p,a = 

= ^{{V? + Cf + 2VlCl){Vs + Cs))p,a = 

= ^iK' + {Cf)p,a)Vi + 2ViPu + P-f{GfGs)p,a- (63) 

Substituting Eq. (62) in Eq. (57), and substituting Eqs. (61), (63) in Eq. (58), we get the 
equations for the functions pM{<l,t), ^^(q, t), £^(q, t). From (55), we get the continuity equation for 
fractional systems [17, 18]: 

-Q-^PM + q^PmVs = 0, (64) 

which may be regarded as the equation of balance of "fractional matter" . This matter can be described 
by fractional systems (30). In addition to pM, the continuity equation (64) includes the density of 
momentum pM^ ■ To obtain the equation for the density of momentum, we multiplied first fractional 
Bogoliubov equation by p°, and use fractional integration over p. Taking advantage of the assumption 



(53) and the boundary condition (51), we get the equation for the components of the vector of density 
of momentum 

-JmVi + —{PmVsVi + Psi) = fi, (65) 



dt' 



where fi = //(q,t) = Z{Fs)p^a- 

Finally, to write down the equation of balance of fractional kinetic energy density, we multiplied 
first fractional Bogoliubov equation by p", and use fractional integrated with respect to p. Taking 
advantage of the assumption (53) and the boundary condition (51), we come to the equation for the 
density of kinetic energy of fractional systems: 



U¥^--'^^)-ik{'' 



PmV^ 



+ E 



+ PsiVi + Qs 



fsVs, 



(66) 



where 



Qs — -rri^i ^s)p,a- 



Equations (64), (65), and (66) are the hydrodynamic equations for fractional systems. Obviously, the 
set of this five equations is not closed. If we have tTsI = and Qs = 0, then these equations are 



d d dP 

^PMVl + ^iUMV.V,) = -— + /„ 



l(5^"''^ + ^)+|fl''- 



If 



PmV^ 



+ E + P 



fsVs, 



(67) 
(68) 
(69) 



and the set of equations is closed. 

7 Conclusion 

In this paper the hydrodynamic equations for fractional systems are derived. In order to derive 
these equations, we use the first Bogoliubov equation for fractional systems [18]. Then we define the 
fractional generalization of the average values and the reduced distribution functions. The Enskog 
equation for fractional systems is considered. Hydrodynamic equations (64), (65), and (66) can be 
considered as equations in the fractional space [17, 18] or for systems with non-Gaussian statistics 
[17, 18]. 

Dissipative and non-Hamiltonian systems can have stationary states of the Hamiltonian systems 
[39]. Classical dissipative systems can have canonical Gibbs distribution as solutions of Liouville 
equations for the dissipative systems [40, 41, 14]. Using the methods [40, 14] it is easy to find 
solutions for the Bogoliubov equations for fractional systems. Suggested Bogoliubov equation allows 
to formulate the dynamics for fractional generalization of quantum dissipative systems by methods 
suggested in [36, 37, 38]. 
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